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ABSTRACT 

Number counts of galaxies are re-analyzed using a semi-analytic model (SAM) of galaxy formation based on 
the hierarchical clustering scenario. Faint galaxies in the Subaru Deep Field (SDF, near-infrared J and K') and 
the Hubble Deep Field (HDF, ultraviolet/optical U , B, V, and /) are compared with our model galaxies. We have 
determined the astrophysical parameters in the SAM that reproduce observations of nearby galaxies, and used 
them to predict the number counts and redshifts of faint galaxies for three cosmological models, the standard 
cold dark matter (CDM) universe, a low-density flat universe with nonzero cosmological constant, and a low- 
density open universe with zero cosmological constant. The novelty of our SAM analysis is the inclusion of 
selection effects arising from the cosmological dimming of surface brightness of high-redshift galaxies, and from 
the absorption of visible light by internal dust and intergalactic H I clouds. As was found in our previous work, 
in which the ultraviolet/optical HDF galaxies were compared with our model galaxies, we find that our SAM 
reproduces counts of near-infrared SDF galaxies in a low-density universe either with or without a cosmological 
constant, and that the standard CDM universe is not preferred, as suggested by other recent studies. Moreover, 
we find that simple prescriptions for (1) the timescale of star formation being proportional to the dynamical time 
scale of the formation of galactic disks, (2) the size of galactic disks being rotationally supported with the same 
specific angular momentum as that of surrounding dark halo, and (3) the dust optical depth being proportional to 
the metallicity of cold gas, cannot completely explain all of observed data. Improved prescriptions incorporating 
mild redshift-dependence for those are suggested from our SAM analysis. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - large-scale structure of 
universe 



1. INTRODUCTION 

It is well known that the number of faint galaxies in a given 
area of sky can constrain cosmological parameters, because it 
depends on the geometry of the Universe (e.g., Peebles 1993). 
Many efforts have been devoted to this subject using traditional 
galaxy evolution models assuming monolithic collapse, such as 
the wind model for elliptical galaxies and the infall model for 
spiral galaxies (e.g., Yoshii & Takahara 1988). These models 
are in fact able to reproduce many of the observed properties 
of nearby galaxies, and provide a useful theoretical tool for un- 
derstanding their evolution (Arimoto & Yoshii 1986, 1987; Ari- 
moto, Yoshii & Takahara 1991). 

In the analyses of galaxy counts from traditional approaches, 
it has been found that the Einstein-de Sitter (EdS) universe, 
a representation of the standard cold dark matter (CDM) uni- 
verse, is not reconcilable with the observed high counts of faint 
galaxies, and that a low-density universe is preferred (Yoshii 
& Takahara 1988; Yoshii & Peterson 1991, 1995; Yoshii 1993). 
Recently, Totani & Yoshii (2000) and Totani et al. (2001a) com- 
pared their predictions against the observed number counts to 
the faint limits in the Hubble Deep Field (HDF; Williams et al. 
1996) and in the Subaru Deep Field (SDF; Maihara et al. 2001), 
taking into account various selection effects and allowing for 
the possibility of number evolution of galaxies in a phenomeno- 
logical way. Note that the SDF counts are now the deepest near- 
infrared ones with the Sa limiting magnitude of K = 23.5 in to- 



tal magnitude. They confirmed that the EdS universe cannot 
reproduce the observed high counts. However, in their best-fit 
models, the merger rates of HDF and SDF galaxies are a little 
different. A mild merger rate is needed to reproduce the counts 
in the HDF, while a negligible rate was necessary for the SDF. 
It should be noted that the photometric passbands for the two 
applications are different; ultraviolet/optical for the HDF and 
near-infrared for the SDF. They suggested that the difference of 
the merger rate might be originated by morphology-dependent 
number evolution because late-type galaxies are mainly seen 
in shorter wavelength such as B-band and early-type galaxies 
are seen in longer wavelength such as /T-band. In any case, it 
should be explained why the merger rate depends on the ob- 
served wavelength in order to obtain a better understanding of 
the galaxy formation process. 

On the other hand, in the study of formation of large-scale 
structure in the universe, both theory and observation suggest 
that gravitationally bound objects, such as galaxy clusters, are 
formed through continuous mergers of dark halos with an ini- 
tial density fluctuation spectrum predicted by the CDM mod- 
els. Based on this scenario of hierarchical clustering, the so- 
called semi-analytic models (SAMs) of galaxy formation have 
been developed by a number of authors (Kauffmann, White & 
Guiderdoni 1993; Cole et al. 1994, 2000; Somervifle & Pri- 
mack 1999; Nagashima et al. 2001, hereafter NTGY). SAMs 
successfully reproduced a variety of observed features of local 
galaxies, such as their luminosity function, color distribution. 
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and so on. 

Faint galaxy number counts have also been analyzed us- 
ing SAMs (Cole et al. 1994; Kauffmann, Guiderdoni & White 
1994; Heyl et al. 1995; Baugh, Cole & Frenk 1996). These 
studies showed that predicted number counts in the EdS uni- 
verse agree with the observed counts. Their results, however, 
contradict analyses carried out with traditional galaxy evolution 
models. Recently, this contradiction was resolved by NTGY, in 
which their SAM is compared with the galaxy counts in the 
HDF. They found, by matching properties of model galaxies 
with observation especially in local luminosity functions and 
cold gas mass fraction, that normalization of model parame- 
ters, related to combinations of physical processes such as star 
formation (SF) and supernova (SN) feedback, is very impor- 
tant, and that accounting for selection effects caused by cosmo- 
logical dimming of surface brightness and absorption of emit- 
ted light by internal dust and intergalactic H I clouds is crucial 
in the analysis of galaxy counts, as shown by Totani & Yoshii 
(2000). It should be noted that recent analysis by Lanzetta et 
al. (2001) also clarify the importance of the selection effects 
caused by the cosmological dimming of surface brightness in 
the observational point of view. They introduced the star for- 
mation rate intensity distribution function, which was derived 
from the ultraviolet luminosity density for the HDF galaxies, at 
several redshifts and found that at high redshift significant frac- 
tion of ultraviolet luminosity density must be missed due to the 
cosmological dimming of surface brightness. 

The purpose of this paper is to examine whether our SAM 
can simultaneously reproduce both the ultraviolet/optical and 
near-infrared galaxy counts in the HDF and in the SDF. Because 
the luminosity of galaxies in different passbands reflects the 
influence of different stellar populations, the subject of multi- 
band number counts provides a strong constraint on galaxy for- 
mation. In this paper, using selection criteria for SDF galaxies 
based on Totani et al. (2001a), we compare our SAM prediction 
of galaxy counts with the observed counts in the SDF. 

This paper is outlined as follows. In §2 we briefly describe 
our SAM, which is almost the same as our previous models 
(NTGY). In §3 we constrain the astrophysical parameters in 
our SAM analysis using local observations. In §4 we compare 
theoretical number counts of faint galaxies with the HDF and 
SDF data, and in §5 we discuss the range of uncertainties in 
our calculations of galaxy number counts. In §6 we provide a 
summary and discussion. 

2. MODEL 

The galaxy formation scenario that we model is as follows. 
In the CDM universe, dark matter halos cluster gravitationally, 
and merge in a manner that depends on the adopted power spec- 
trum of the initial density fluctuations. In each of the merged 
dark halos, radiative gas cooling, star formation, and gas re- 
heating by supernovae occur. The cooled dense gas and stars 
constitute galaxies. These galaxies sometimes merge together 
in a common dark halo, and more massive galaxies form. Re- 
peating these processes, galaxies form and evolve to the present 
epoch. 

The SAM analysis we perform obtains essentially the same 
results of our previous SAM analyses, with minor differences 
in a number of details. In this section we briefly describe our 
model. Details of the model we employ are described in NTGY, 
except for a few differences, which are explicitly mentioned be- 
low. 



2.1. Scheme of Galaxy Formation 

The merging histories of dark halos are realized by a 
Monte Carlo method proposed by Somerville & Kolatt (1999), 
based on the extended Press-Schechter formaUsm (Bond et 
al. 1991; Bower 1991; Lacey & Cole 1993). This formal- 
ism is an extension of the Press-Schechter formalism (Press & 
Schechter 1974), which gives us the mass function of dark ha- 
los, n{M)dM, to estimate the mass function of progenitor halos 
with mass M\ at a redshift zo + Az of a single dark halo with 
mass Mo collapsing at a redshift zo, n(Mi;zo + Az|Mo;zo)t^Mi. 
Dark halos with circular velocity Vcirc > 40km s"' are regarded 
as isolated halos, and otherwise regarded as diffuse accreted 
matter. The mean density in dark halos is assumed to be pro- 
portional to the cosmic mean density at the collapsing epoch 
using a spherically symmetric collapse model (Tomita 1969; 
Gunn & Gott 1972). Each collapsing dark halo contains bary- 
onic matter with a mass fraction f2b/^^o> where fio and fib are 
the parameters of total and baryon mass densities relative to the 
critical cosmic mass density. We adopt a value of fib = 0.015/?"^ 
given by Suzuki, Yoshii & Beers (2000) in which they estimated 
the value from the primordial lithum abundance, where h is the 
Hubble parameter, h = HQ/ 100km s"^Mpc"^ . Note that a recent 
measurement of the anisotropy of the cosmic microwave back- 
ground by the BOOMERANG project suggests a slight higer 
value, fib ^ 0.02/!-2 (Netterfield et al. 2001). The effect of 
changing fib has already been investigated by Cole et al. (2000) 
and they showed that this mainly affects the value of the in- 
visible stellar mass fraction such as brown dwarfs parametrized 
by T (see the next subsection). We also checked whether our 
results are changed or not in the case of fib ~ 0.04 and found 
that this does not affect them. The baryonic matter consists of 
diffuse hot gas, dense cold gas, and stars. 

When a halo collapses, the hot gas is shock-heated to the 
virial temperature of the halo with an isothermal density pro- 
file. A part of the hot gas cools and accretes to disk of a galaxy 
until the subsequent collapse of the dark halos. The amount 
of the cold gas involved is calculated by using metalUcity- 
dependent cooling functions provided by Sutherland & Dopita 
(1993). The difference of cooling rates between the primor- 
dial and metal-polluted gases is prominent at J ~ lO^K due to 
Une-cooUng of metals. The cooUng is, however, very efficient 
in dark halos with a virial temperature of T ^ lO^K even in 
the case of the primordial gas, so the metallicity dependence of 
cooling rate only slightly affects our results. In order to avoid 
the formation of unphysicaUy large galaxies, the cooUng pro- 
cess is applied only to halos with Vcirc <500 km s"^ in the stan- 
dard CDM and 400 km s"' in low-density universes according 
to some previous SAMs. This handling would be needed be- 
cause the simple isothermal distribution forms so-called "mon- 
ster galaxies" due to the too efficient cooling at the center of 
halos. While this problem will probably solved by adopting 
another isothermal distribution with central core (Cole et al. 
2000), we take the above simple approach because this does 
not cause any serious problems in estimating galaxy counts. 

Stars in the disk are formed from the cold gas. The SF 
rate (SFR) is given by the cold gas mass Mcou and a SF 
timescale r* as = McoIcj/t,. We use two SF models. One 
model adopts constant star formation (CSF), in which r* is a 
constant against redshift; the other adopts dynamical star for- 
mation (DSF), in which is proportional to the dynamical 
timescale of the halo, which allows for the possibility that the 
SF efficiency is variable with redshift. We then express the SF 
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timescale as 
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(CSF), 



(l+z)'^(DSF), 



(1) 



where T°,a* and cr are free parameters. The former two pa- 
rameters are chosen so as to match the fraction of cold gas 
mass with the observed fraction (see next section). While, of 
course, some fraction of cold gas might not be seen, we adopt 
the gas mass fraction to constrain the SF-related parameters ac- 
cording to Cole et al. (2000). If much of this cold gas would 
be hidden, should be longer than our value. Note that in 
our previous paper (NTGY) we used DSF only in the case of 
ij = 0. In Figure I we show the timescales for CSF (dot-dashed 
line) and DSF with cr = (solid line) and (7=1 (dashed line) in 
a A-dominated flat universe. Because objects which collapsed 
at higher redshift have higher density, the dynamical timescale 
Tdyn is shorter at higher redshift. Thus, the DSF scenario with 
smaller a converts the cold gas into stars more rapidly than the 
CSF. In this paper we examine not only the CSF and cr = DSF 
models but also the DSF model with a = 1 as an intermediate 
case. 
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Fig. 1. — Star formation timescales r, as a function of redshift z in a A- 
dominated flat CDM model. The timescale is normalized by that evaluated at 
present. The dot-dashed line denotes the constant SF model (CSF). The solid 
line denotes the dynamical star formation model (DSF) with (t = for which 
the timescale is proportional to the dynamical timescale of disk. The dashed 
line denotes the DSF with it = 1 as an intermediate case between CSF and DSF 
with (T = 0. 

Massive stars explode as Type II supemovae (SNe) and heat 
up the surrounding cold gas. This SN feedback reheats the cold 
gas at a rate of Mreheat = ^coid/7reheat> where the timescale of 
reheating is given by 



"^reheat ~ 



Ki. 
Vho, 



r*. 



(2) 



The free parameters Vh„t and cthot are determined by matching 
the local luminosity function of galaxies with observations. 

With the parameters M* and Mreheat thus determined, we ob- 
tain the masses of hot gas, cold gas, and disk stars as a func- 
tion of time during the evolution of galaxies. Chemical enrich- 
ment is also taken into account, adopting heavy-element yield 
of y = 0.038 = 2Z0, but changing this value of y has a minimal 
effect on the results described below. 



When two or more progenitor halos have merged, the newly 
formed larger halo should contain at least two or more galaxies 
which had originally resided in the individual progenitor halos. 
By definition, we identify the central galaxy in the new com- 
mon halo with the central galaxy contained in the most mas- 
sive one of the progenitor halos. Other galaxies are regarded 
as satellite galaxies. These satellites merge by either dynami- 
cal friction or random collision. The timescale of merging by 
dynamical friction tmc is given by Binney & Tremaine (1987), 
which is estimated from the baryonic mass of satellite. When 
the time elapsed after merging of progenitor halos exceeds tmc, 
a satellite galaxy in the common halo is accreted to the cen- 
tral galaxy. On the other hand, the mean free timescale of ran- 
dom collision Tcou is given by Makino & Hut (1997). With a 
probability Af / Tcou, where Af is the time step corresponding 
to the redshift interval Az of merger tree of dark halos, a satel- 
lite galaxy merges with another satelUte picked out randomly. 

Consider the case when two galaxies of masses m\ and 
m2 (> ffJi ) merge together. If the mass ratio f = m\lmi\<s, larger 
than a certain critical value of /bulge, we assume that a starburst 
occurs and that aU of the cold gas turns into stars and hot gas, 
which fills the resulting halo, and all of the stars populate the 
bulge of a new galaxy. On the other hand, if / < /bulge, no 
starburst occurs, and a smaller galaxy is simply absorbed into 
the disk of a larger galaxy. This division into major and mi- 
nor modes is only for simpUcity and we adopt /bulge =0.1 as a 
standard value. As shown in Figure 13, the dependence of to- 
tal counts on /bulge is negligible (see §5.4). Although there is a 
model to treat continuous burst activity, say, as a function of the 
mass ratio / (Somerville, Primack & Faber 2001), we adopt the 
above simple division because of the negligible denpendence 

on /bulge- 

2.2. Photometric Properties of Galaxies 

The above processes are repeated until the output redshift 
and then the SF history of each galaxy is obtained. For the pur- 
pose of comparison with observation, we use a stellar popula- 
tion synthesis approach, from which the luminosities and colors 
of model galaxies are calculated. Given the SFR as a function 
of time or redshift, the absolute luminosity and colors of indi- 
vidual galaxies are calculated using a population synthesis code 
by Kodama & Arimoto (1997). The stellar metalUcity grids in 
the code cover a range from Z* =0.0001 to 0.05. Note that we 
now define the metallicity as mass fraction of metals, for ex- 
ample, the solar metallicity is 0.019. The initial stellar mass 
function (IMF) that we adopt is the power-law IMF of Salpeter 
form, with lower and upper mass limits of O.IM© and 6OM0, 
respectively. Then, following Cole et al. (1994), we introduce 
a parameter defined as T = (Mjum + Mbd) /Mum, where Mium is 
the total mass of luminous stars with m > Q.IMq and Mbd is 
that of invisible brown dwarfs. 

Our model of estimating the optical depth of internal dust 
has been improved over that adopted by NTGY. We take the 
usual assumption that the abundance of dust is proportional to 
the metallicity of cold gas, and that the optical depth is pro- 
portional to the colunm density of metallicity. Then the optical 
depth r is given by 

Mcold-Zcold , 

rex 2 (1+z) ' (3) 

where r^ is the effective radius of the galactic disk, described in 
the next subsection. The proportionaUty constant is determined 
by matching the extinction. Ay ~ 0.2, for galaxies typical of the 
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Milky Way, assuming the slab dust model according to our pre- 
vious paper (NTGY) and Somerville & Primack (1999). We ad- 
just the high redshift properties of dust-to-metal ratio and dust 
dumpiness by changing 7. As a standard model, we simply 
adopt 7 = 0. It will be also investigated in §5.3 how the galaxy 
counts are affected by 7. 

For some classes of galaxies, the screen dust model has been 
suggested to be a good approximation, rather than the slab dust 
model. For example, some galaxies such as nearby starburst 
galaxies and hyper extremely red objects (HEROs) are very red 
and therefore their dust distribution should be approximated by 
the screen model (Totani et al. 2001b). Besides recent analysis 
of source counts in far-infrared and submilimeter wavelengths 
also favors the screen dust model (Totani & Takeuchi 2002). 
In this paper, however, we adopt the slab model as a standard 
one because main contribution in UV/optical and near-infrared 
passbands is considered to be galaxies which are not heavily 
extincted under modest star formation activity. Of course, even 
for these normal galaxies the screen dust model might be a good 
approximation. Therefore we will check the screen dust model 
in §5.3. While we check only these two models, the uncertainty 
cased by the dust model would be clarified by checking these 
models and by varying 7. 

Emitted hght from distant galaxies is absorbed by Lyman 
lines and Lyman continuum in intervening intergalactic H I 
clouds. We used an optical depth calculated by Yoshii & Pe- 
terson (1994) to account for this. 

We classify galaxies into different morphological types ac- 
cording to the fi-band bulge-to-disk luminosity ratio B/D. In 
this paper, following Simien & de Vaucouleurs (1986), galaxies 
with B/D > 1 .52, 0.68 < B/D < 1 .52, and B/D < 0.68 are clas- 
sified as ellipticals, SOs, and spirals, respectively. Kauffmann, 
White & Guiderdoni (1993) and Baugh, Cole & Frenk (1996) 
showed that this method of type classification reproduces the 
observed type mix well. 

We then assess whether the surface brightnesses of model 
galaxies are above the detection thresholds of the SDF and HDF 
observations. The selection effects in our SAM analysis are 
evaluated as follows. Using the intrinsic size of model galax- 
ies, described in the next subsection, their surface brightness 
profiles in the observers frame are given by a convolution of 
the luminosity profile with a Gaussian point-spread function. 
Here we assume an exponential profile for spiral galaxies, and 
a de Vaucouleurs profile for elUptical and lenticular galaxies. 
We note that model galaxies with surface brightness higher than 
the threshold 5th, and isophotal areas larger than the minimum 
area Ath are actually detected as galaxies (Yoshii 1993). Fol- 
lowing Totani & Yoshii (2000) and Totani et al. (2001a), we 
adopt 5th = 25.59 mag arcsec"^ in K' , 5th = 24.10 mag arcsec"^ 
in J, and Ath =0.24 arcsec^ for the SDF galaxies, and 5th = 27.5 
mag arcsec"^ in Veoe, 5th = 27.0 mag arcsec"^ in /8i4,f/3oo and 
B450, and Ath =0.04 arcsec^ for the HDF galaxies. 

For the SDF galaxies, we also take into account the detection 
probability, that is, completeness. Noise and statistical fluctua- 
tions in the data prevent complete source detections at any spec- 
ified limits of surface brightness and size. Totani et al. (2001a) 
found that the dispersion of observed isophotal area around its 
true value can be fitted by 



aA{ni,dob) = c(Ai -A2)'/^c/oi 



(4) 



the isophotal level 0.8 and 1 .2 times brighter than 5th, respec- 
tively. Assuming a Gaussian distribution of observed isophotal 
area with this dispersion, we obtain the probabiUty that the ob- 
served isophotal area is larger than the threshold value Ath. Ad- 
ditional details of these selection effects are described in Totani 
& Yoshii (2000) and Totani et al. (2001a). 

2.3. Galaxy Size 

We assume that the size of spiral galaxies is determined by a 
radius at which the gas is supported by rotation, under the con- 
servation of specific angular momentum of hot gas that cools 
and contracts. We also assume that the initial specific angular 
momentum of the gas is the same as that of the host dark halo. 
Acquisition of the angular momentum of dark halos is deter- 
mined by the tidal torques in the initial density fluctuation field 
(White 1984; Catelan & Theuns 1996a,b; Nagashima & Gouda 
1997). The distribution of the dimensionless spin parameter Ah 
is well approximated by a log-normal distribution (Mo, Mao & 
White 1998), 



1 



liTcrx 



exp 



(InAH-lnA)^ 



6/ In Ah, (5) 



where m is the total magnitude of an object, dob is the FWHM 
size, and Ai and A2 are the isophotal areas corresponding to 



where A is the mean value of spin parameter and a\ is its loga- 
rithmic variance. We adopt A = 0.05 and a\ = 0.5. (Note that in 
our previous paper we adopted Ah = 0.05 for all spirals.) If the 
specific angular momentum is conserved, the effective radius 
of a presently observed galaxy at z = is related to the initial 
radius r, of the progenitor gas sphere via = (1.68/a/2)Ah?"; 
(Fall 1979; Fall & Efstathiou 1980; Fall 1983). The initial ra- 
dius r, is set to be the smaller one between the virial radius of 
the host halo and the cooling radius. A disk of a galaxy grows 
due to cooling and accretion of hot gas from more distant enve- 
lope of its host halo. In our model, when the estimated radius 
by the above equation becomes larger than that in the previous 
step, the radius grows to the new larger value in the next step. 

Size estimation of high-redshift spiral galaxies, however, car- 
ries uncertainties because of the large dispersion in their ob- 
served size distribution. Allowing for the possibility that the 
conservation of angular momentum is not complete, we gener- 
alize this size estimation by introducing a simple redshift de- 
pendence, 

re = ^XHri{l+zr, (6) 

where p is a free parameter; we simply use p = as a refer- 
ence value in this paper. Results for other values of p will be 
given in §5. The effect of changing p emerges in the selection 
effects due to the cosmological dimming of surface brightness 
and in the dust extinction because the dust column density also 
changes with galaxy size. 

Figure 2 shows the effective disk radii of spiral galaxies 
at z = as a function of their luminosity for the models of 
SC (thick sohd line), OC (dashed Une), and LC (dot-dashed 
line). The adopted parameters in these models are tabulated 
in Table 1. All the models well reproduce the observed disk 
size-magnitude relation (thin sohd line) compiled by Totani & 
Yoshii (2000) based on the data taken from Impey et al. (1996). 
The redshift dependence of such relation through changing p in 
equation 6 will be discussed in §5. 

The sizes of early-type galaxies, which likely form from 
galaxy mergers, are determined by the virial radius of the bary- 
onic component. When a major merger of galaxies occurs, the 
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merged system consisting of stars and cold gas initially virial- 
izes to satisfy 



2 ^ G(M^ +Mcoid) 

circ 



(7) 



In this paper we take into account the dynamical response to 
the "shallowing" of the gravitational potential due to the mass 
loss caused by the SN feedback. Then, after the initial virial- 
ization, the system expands while losing the cold gas, with the 
adiabatic invariant (M, +Mcoid)'' kept constant (e.g., Yoshii & 
Arimoto 1987). Then the final radius of the system, consisting 
only of stars after the loss of the cold gas, is given by 



'"/= 



M, 



(8) 



J 



where M*,, and M* j are the total masses of stars in the sys- 
tem before and after the mass loss, respectively. The effect 
of the dynamical response is most prominent for dwarf galax- 
ies of lower circular velocity, and can explain the properties 
of local dwarf ellipticals, which will be presented in a sepa- 
rate paper. However, this effect has only a minor impact on 
our predicted galaxy number counts, because luminous galax- 
ies of larger circular velocity are the main contributors to galaxy 
number counts at the faint limits under consideration. In order 
to match with the observed effective radius, we here introduce 
a scahng parameter, /b, such that = f^rf. The adopted values 
of /b for SC, OC, and LC are tabulated in Table 1 . 
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Fig. 2. — Galaxy sizes. The thick solid line connecting open circles, 
the dashed line connecting crosses, and the dot-dashed line connecting open 
squares show theoretical results for effective radii of late-type galaxies in the 
SC, OC, and LC models, respectively. We show Icr scatter in the LC model 
by error bars. Similar scatter in other models are not shown. The thin line 
shows the observed best-fit relation given by Totani & Yoshii (2000) for spiral 
galaxies, based on the data taken from Impey et al. (1996). 



3. PARAMETER SETTINGS 

We normalize the model parameters so as to agree with vari- 
ous local observations. As shown in NTGY, this normalization 
procedure is an essential ingredient in the analysis of galaxy 
counts. 

The cosmological parameters (f2o> ^a, h, fg) adopted in this 
paper are tabulated in Table 1. For all the models (SC, OC, 
and LC), the baryon density parameter fib = 0.015/;"^ is used 
in common. For OC and LC the value of erg = 1 is determined 
from observed cluster abundances (Eke, Cole & Frenk 1996). 

The astrophysical parameters (Vhot, ahot, t^, a*, fi,, /bulge. 
T) are constrained from local observations as discussed below. 



The adopted values are almost the same as in our previous pa- 
per. Slight differences have resulted from modifications to the 
model, such as the improved scheme of dust extinction. In this 
paper we only show the result of CSF (equation 1) as a stan- 
dard, and other parameters used in common, including 7 = 
(equation 3), p = (equation 6), and (A,(Ta)=(0.05, 0.5) (equa- 
tion 5). 
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Fig. 3 . — Local luminosity functions in the (a) B band and (b) K band. The 
solid, dashed, and dot-dashed lines indicate the SC, OC, and LC models, re- 
spectively. Symbols with error bars in (a) indicate the observational data from 
APM (Loveday et al. 1992, filled squares), ESP (Zucca et al. 1997, open tri- 
angles), DurhamAJKST (Ratcliffe et al. 1998, open diamonds), 2dF (Folkes et 
al. 1999, open circles), and SDSS (Blanton et al. 2000, crosses). Note that the 
SDSS limiinosity function shown here is that with the same detection limit as 
employed in the 2dF survey. Symbols in (b) indicate the data from Mobasher 
et al. (1993, open squares), Gardner et al. (1997, open stars), and 2MASS 
(Cole et al. 2000, filled triangles). 

First, the SN feedback-related parameters (Vhot) cthot) and the 
mass fraction T of invisible stars are almost uniquely deter- 
mined if their values are chosen so as to reproduce the local 
luminosity function. Figure 3 shows theoretical results for the 
models of SC (thick solid Une), OC (dashed Une), and LC (dot- 
dashed line), which are tabulated in Table 1 . Symbols with er- 
ror bars indicate observational results from the B-band redshift 
surveys, such as APM (Loveday et al. 1992), ESP (Zucca et 
al. 1997), Durham/UKST (Ratchffe et al. 1998), 2dF (Folkes 
et al. 1999) and SDSS (Blanton et al. 2001), and from the K- 
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band redshift surveys (Mobasher et al. 1993; Gardner et al. 
1997; 2MASS, Cole et al. 2001). As done in NTGY, we regard 
the APM result as a standard one. This gives somewhat less 
normalization of luminosity function compared with other sur- 
veys', probably because of its worse completeness compared 
to the other recent surveys. Nevertheless, in this paper, we de- 
termine the SN feedback-related parameters in the above way 
because the dependence of these parameters on galaxy counts 
has been already investigated by NTGY and then because it has 
been found that the shape of counts is not significantly affected. 
We checked that the shape of the counts in the /T'-band are also 
not affected by changing the luminosity function significantly. 

Next, the SFR-related parameters (t^, a,) are determined by 
using the cold gas mass fraction in late-type galaxies. The gas 
fraction depends on both the SN feedback-related parameters 
and on the SFR-related ones. The former parameters determine 
the gas fraction expelled from galaxies and the latter ones the 
gas fraction that is converted into stars. Therefore, in advance 
of determining the SFR-related parameters, the SN feedback- 
related parameters must be determined by matching the local 
luminosity function. 

Figure 4 shows the ratio of cold gas mass relative to B-band 
luminosity of spiral galaxies as a function of their luminosity. 
Theoretical results are shown for the models of SC (thick sohd 
line), OC (dashed line), and LC (dot-dashed line). We here as- 
sume that 75% of the cold gas in these models is comprised of 
hydrogen, i.e., Mhi = 0.75Mcoid. H I data, taken from Hucht- 
meier & Richter (1988), are shown by open squares with er- 
ror bars. Since their data do not include the fraction of H2 
molecules, the observational result should be regarded as pro- 
viding a lower Umit of the cold gas mass fraction. 
T ' r 
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Fig. 4. — Cold gas mass relative to B-band luminosity of spiral galaxies. 
The solid, dashed, and dot-dashed curves show the SC, OC, and LC models, 
respectively. The open squares indicate the observational data for atomic neu- 
tral hydrogen taken from Huchtmeier & Richter ( 1 988). In the models the cold 
gas consists of all species of elements, therefore its mass is multiplied by 0.75, 
i.e.. A/hi = 0.75Mcoid, which corresponds to the fraction of hydrogen. Because 
the observational data denote only atomic hydrogen, they should be regarded 
as lower limits of the ratio. 

4. RESULTS 

4.1. Galaxy Number Counts 

Figure 5 shows the galaxy number counts in the Subaru K'- 
band as a function of apparent isophotal K' magnitude for the 
models of SC, OC, and LC. Thick lines show theoretical predic- 
tions based on the observational conditions, including the selec- 
tion effects from the cosmological dinnming of surface bright- 



ness and completeness for galaxies in the SDF Open circles 
indicate the observed SDF raw counts in isophotal magnitude. 
The number of galaxies in the SC model is too small to ex- 
plain the observed raw counts (open circles with error bars). 
On the other hand, the OC and LC models explain the SDF 
counts equally well. For reference, the data from other surveys 
are also shown by the specified symbols, after applying a trans- 
formation of K' = K+0.1. 

Thin lines in this figure are theoretical predictions without 
the selection effects. It is evident that the effects become impor- 
tant at K' > 20, and make the count slope turn over at K' ~ 24, 
where the observed raw counts reach a maximum. The pre- 
dicted counts with and without the selection effects differ by a 
factor of three at K' ^ 24, which is much greater than the ob- 
servational error. Thus, incorporation of the selection effects is 
essential in the SAM analysis of the galaxy number counts. 

We check for consistency with the HDF counts in the UBVI 
bands in Figure 6. The types of lines are the same as in Fig- 
ure 5. We see that the OC and LC models, which reproduce 
the near-infrared SDF counts, give a better agreement with 
the UV/optical HDF counts when compared with our previous 
models in NTGY. On the other hand, the SC model significantly 
falls below those observed in all the UBVIK' bands. Thus we 
conclude that our SAMs well reproduce the multiband galaxy 
number counts in a low-density universes with or without a cos- 
mological constant. 

4.2. Isophotal Area-Magnitude Relation 

As stressed in the previous section and in NTGY, the selec- 
tion effects from the cosmological dimming of surface bright- 
ness of galaxies cannot be ignored in the SAM analysis of 
galaxy number counts. This indicates that the size of high- 
redshift galaxies must be modeled properly. Figure 7 plots the 
SDF galaxies in the isophotal area-magnitude diagram. The 
data plotted are those for the SDF galaxies that are detected in 
both the K' and J bands. Only the LC model is shown, with var- 
ious values of p = (solid line), 0.5 (dashed line), and 1 (dot- 
dashed line), because other parameters involving the adopted 
cosmology and dust extinction hardly affect the result. The 
isophotal area becomes smaller for fainter apparent magnitude; 
at faint Umits it reaches the minimum for detection. All pre- 
dictions with p = - 1 give a convergent result at faint limit of 
K' ^ 24. This convergence occurs because galaxies with larger 
areas at K' ~ 24 have surface brightnesses below the threshold 
and remain undetected. We find from this figure that our SAM 
galaxies well reproduce the observed area-magnitude relation, 
and are consistent with the SDF galaxies, only when the se- 
lection bias against faint galaxies with high redshift and/or low 
surface brightness is taken into account in the analysis. 

4.3. Redshift Distributions 

The left panels in Figure 8 show the redshift distributions for 
the HDF galaxies in the I^u band. The sohd and dashed Unes 
are theoretical predictions for the LC model with and without 
the selection effects, respectively. The histograms are the ob- 
served photometric redshift distributions given by Furusawa et 
al. (2000). As sown in NTGY, the LC model well reproduces 
the observed redshift distribution as a function of apparent mag- 
nitude. 

The right panels of Figure 8 show the /T'-band predictions for 
the LC model. Compared with Totani et al. (2001a), in which 
the traditional models of galaxy evolution are used, the peak 
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location for our redshift distribution moves to lower redshift. 
Therefore, redshifts of the SDF galaxies, if measured by either a 
spectroscopic or photometric method, give an important insight 
into physical processes of galaxy formation and evolution. 

5. MODEL UNCERTAINTIES AT HIGH REDSHIFT 

In this section we discuss uncertainties in predicting the num- 
ber counts of galaxies. Such uncertainties arise mainly from the 
SF timescale (§5.1), the galaxy size (§5.2), and the dust extinc- 
tion (§5.3). As for the SN feedback-related parameters, we find 
that the dependence of our result on their values has a similar 
tendency to NTGY, so that our conclusion of preferring a low- 
density universe is not altered by uncertainty associated with 
these parameters. Finally we will discuss the morphological 
counts (§5.4). 

5.1. Star Formation Timescale 

We evaluate the effects of changing the SF timescale from 
CSF to DSF (fj = or 1) in the ACDM model (see equation 1). 
Note that CSF in the EdS universe is equivalent to DSF with 
(7 = 1.5, thus the case of DSF with cr = 1 is in between CSF and 
DSF with a = 0, as shown in Figure 1. Hereafter, the two cases 
of DSF with (7 = and 1 in the ACDM model are referred to as 
LDO and LDl, respectively. 

Figure 9 shows the predicted number counts of galaxies in 
the Subaru K' and HST B450 bands for the models of LC (solid 
line), LDl (dashed line) and LDO (dot-dashed line). The pre- 
dicted counts from LDl and LDO are higher than those by LC; 
this tendency is much more prominent in the longer wavelength 
passband. This wavelength-dependent tendency is explained 
as follows. In general, the near-infrared luminosity in the K' 
band mainly reflects the total mass of long-hved stars, while 
the UV/blue luminosity reflects the instantaneous SFR. While 
the total mass of cold gas available for star formation is limited 
by the radiative cooling of hot gas, the conversion timescale 
from cold gas to stars is determined by the SF timescale. In 
the case of DSF, for which the SF timescale is shorter than in 
CSF, more stars are formed at higher redshift. Consequently, in 
the cases of LDl and LDO, individual galaxies become brighter 
in the K' band, and so the /T'-band counts of observable galax- 
ies above the surface brightness threshold become larger than 
in LC, although their blue luminosity is hardly changed in the 
B450 band. 

The above argument indicates that the SF timescale is best 
constrained by the observed SDF counts in the K' band. Since 
the LC model is found to give a superior agreement with the 
SDF counts, the SF timescale should therefore be nearly con- 
stant against redshift. This result is consistent with previous 
SAM analyses by Kauffmann & Haehnelt (2000), SomerviUe, 
Primack & Faber (2001) and NTGY. 

The optical depth for absorption of stellar light by dust is as- 
sumed to be proportional to the amount of cold gas present. In 
the cases of LDl and LDO, because of the short SF timescale, 
the fraction of cold gas mass at high redshift is smaller, leading 
to smaller optical depth when compared to LC. For this rea- 
son, the B45o-band counts in LDl and LDO, which are sensitive 
to internal dust extinction, show a slight excess at B450 ^ 23. 
However this effect is almost canceled out by the increase of 
stars due to the higher SFR given by the shorter SF timescale, 
as shown in Figure 9b. 

5.2. Galaxy Size 



As mentioned in §4.1, the selection effects from the cosmo- 
logical dimming of surface brightness of galaxies are impor- 
tant. Therefore, the number counts of galaxies are sensitive 
to their size, because the surface brightness is proportional to 
the inverse of radius squared. Figure 10 shows theoretical pre- 
dictions for the LC model with p = Q (solid line), 0.5 (dashed 
line), and 1 (dot-dashed line). Note that the simple model with 
p = leads to smaller disk size at high reshift in proportional 
to the initial radius r,, which is nearly equal to the virial ra- 
dius in the case of halos with galactic mass-scale and becomes 
smaller size according to the spherical collapse model. Thus 
the models with p> lead to larger disks and lower surface 
brightness compared to the simple model with p = 0, and hence 
the smaller number counts of observable galaxies above the sur- 
face brightness threshold, in both the K' and B450 bands, when 
compared to the standard case of p = 0. Since the case of p = 
slightly overpredicts the ^T'-band observed counts at K' ~ 24 
and the case of p ^ 0.5 is preferable, it seems that high-redshift 
galaxies might be more extended in size than expected from the 
cooling radius. 

The size of a galaxy disk also affects dust extinction. The 
optical depth of dust is proportional to the inverse of disk ra- 
dius squared. In the case of p = 1, the number of galaxies 
at £450 ~ 23 increases above the p = prediction, because of 
weaker dust extinction owing to larger disk size. The effects of 
dust extinction are discussed in the next section in detail. 

The physical processes which cause the galaxy disks to be 
more extended at high redshift are yet to be identified. One pos- 
sibility includes the viscous evolution of the star-forming disk 
(Saio & Yoshii 1990), but more investigations will obviously be 
necessary here. 

5.3. Dust Extinction 

The uncertainty in estimation of dust extinction at high red- 
shift comes from a combination of uncertainties in the chemical 
enrichment model, dust-to-metal ratio, dumpiness of dust, and 
so on. Instead of entering into these details, however, we mea- 
sure a robustness of our result by changing the value of 7 in the 
optical depth by dust (equation 3). 

Figure 11 shows theoretical predictions for the LC model 
with 7 = (solid line), 1 (dashed line), and 2 (dot-dashed line), 
assuming the slab dust distribution . The predicted /T'-band 
counts are almost independent of 7 (Figure 11a) because the 
near-infrared passbands are insensitive to dust extinction. On 
the other hand, the predicted B45o-band counts are sensitive to 
dust extinction (Figure lib). In the case of 7 > 1, the "knee" 
in the count curve near B450 ~ 24 mag becomes prominent. It 
might be therefore necessary to adopt a high value of 7 > 1 in 
order to reproduce the observed B45o-band counts better over a 
wide range of B45()-magnitude if the simple assumption of disk 
size in the previous subsection is correct. 

Figure 12 is the same as Figure 11, except for showing the- 
oretical predictions for the screen dust distribution with 7 = 
(dashed line) and 2 (dot-dashed line), in comparison with the 
slab dust distribution with 7 = (solid line). In general, when 
the value of optical depth is fixed in the case of 7 = 0, the screen 
dust model gives stronger extinction than the slab dust model, 
so that the predicted number of bright galaxies at B450 < 23 be- 
comes only slightly lower than that for the slab dust. In the case 
of 7 = 2, the "knee" in the count curve near B450 ~ 24 mag is 
prominent even in the screen dust model. On the other hand, the 
predicted ^'-band counts are not affected at aU, and are there- 
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fore free from the uncertainty in estimation of dust extinction. 

The most effective discriminator of dust extinction is the 
knee of observed galaxy counts at B450 ~ 23. When the dust 
extinction is strong, the knee disappears and the count curve 
deviates below the data. Therefore, the optical depth by dust 
should be small at high redshift, that is, z > 1. We find that the 
knee is sensitive almost exclusively to the optical depth, which 
is determined by the cold gas mass, disk size, metaUicity of cold 
gas, and 7 (equation [3]). 

In order to decrease the optical depth we need to decrease 
the cold gas mass and/or metaUicity, expand the galaxy size, 
or increase the value of 7. If one decreases the cold gas mass 
with a shorter SF timescale, too many galaxies are formed, and 
the number of galaxies in the /T'-band is overpredicted (§5.1). If 
one expands the galaxy size, we suffer from too strong selection 
effects (§5.2). There is an uncertainty in estimation of the cold 
gas metaUicity (Somerville, Primack & Faber 2001), but this 
uncertainty can be practically absorbed in the effect of 7. Thus, 
in the usual SAMs (7 = 0), a simple estimation of optical depth 
might be worse to reproduce the observed galaxy counts at short 
wavelengths rather than models with a high value of 7 > 1 . In 
any case, the shape of galaxy counts at shorter wavelength is 
sensitive to the optical depth. We wiU need more investigations 
on the dust extinction. 

5.4. Morphology-Dependent Number Counts 

Our SAM simultaneously reproduces the multi-band galaxy 
counts, in contrast to the previous works using the traditional 
galaxy evolution models (Totani & Yoshii 2000; Totani et al. 
2001a). The reasons could indeed be quite complicated, but 
one of the possibilities might be that a dominant morphology 
of galaxies is different in different passband. The SAM al- 
lows for morphological transformations during the evolution 
of galaxies, that is, mergers between galaxies of similar mass 
make spheroidals as merger remnants, then they accrete hot gas 
through radiative cooling as galactic disk. 

Figure 13 shows the predicted contributions of early-type 
galaxies (E/SO; dashed lines) in the LC model in the K' and 
/gH bands. Note that the morphology of galaxies is identified 
by the Z?-band bulge-to-disk ratio at their observed redshift in 
the observer frame and that this definition is different from that 
in usual analyses with the traditional models. The total num- 
ber counts are also shown by the soUd Unes. In order to see 
the fractional change of morphologies, we change the value of 
/bulge which divides mergers into major and minor ones accord- 
ing to the mass ratio of merging galaxies (see §2.1). The thick 
lines are for the reference value of /bulge =0.1, and the thin lines 
for /bulge = 0.5 where major merger occurs only when the mass 
of smaller satellite is larger than 0.5 times that of larger one. 
We also plot the observed E/SO counts in the HDF (Abraham et 
al. 1996; Driver et al. 1998; Phillipps et al. 2000) with the total 
HDF counts in the /g^-band. We find that our prediction for 
E/SO counts roughly agree with the observed data and that the 
uncertainty in the E/SO counts caused by that of /bulge is nearly 
a factor of two. 

We note that while the fraction of early-type galaxies in- 
creases with decreasing /bulge, the total counts are hardly 
changed. This indicates that the total number of galaxies is es- 
sentially determined by the merging histories of dark halos and 
by the cooling rate of hot gas, and that the cooling must be very 
efficient to enable frequent formation of disk galaxies at z > 1. 

We also note that our morphological classification is based 



on the B-band bulge-to-disk ratio. It has not yet been con- 
firmed that this 5-band classification coincides with that for the 
K' band, therefore our predictions of type-dependent counts in 
the K' band should be regarded as qualitative estimates. 

6. SUMMARY AND DISCUSSION 

We have calculated the number counts of faint galaxies in the 
framework of SAM for three cosmological models, the stan- 
dard CDM (EdS) universe, a low-density open universe, and a 
low-density flat universe with nonzero A. Our SAM includes 
the selection effects from the cosmological dimming of sur- 
face brightness of galaxies with criteria appropriate for the SDF 
and HDF, and also includes some modifications to our previous 
analysis (NTGY), such as the optical depth estimation of dust 
within a galaxy. In this paper we have shown that our SAM is 
fully consistent with that of the previous version, and can ex- 
plain the observed multi-band galaxy counts from the UV to the 
near-infrared. 

Comparison of theoretical predictions with the observed 
number counts of SDF and HDF galaxies, as well as with other 
ground-based observations, indicates that the standard CDM is 
ruled out, and a A-dominated flat universe and a low-density 
open universe are favored. This result is consistent with that 
from HDF galaxies (NTGY), but is in sharp contrast with pre- 
vious SAM analyses by other authors, where many of the con- 
ceivable selection effects in the faint observations have been 
ignored. 

Some uncertainties in our SAM have been discussed. These 
arise from a lack of knowledge on the galaxy formation pro- 
cess, and also from an insufficient survey of the physical prop- 
erties of high-redshift galaxies. We especially focused on the 
uncertainties in redshift dependence of SF timescale, galaxy 
size, and dust extinction. We found that dust extinction hardly 
affects galaxy counts in the Subaru K' band, but does signif- 
icantly affect those in the HST B450 band. Thus, the ^T'-band 
galaxy counts are robust against the uncertainty of dust extinc- 
tion. Two other factors affect galaxy counts even in the K'- 
band. If the SF timescale at high redshift is shorter than one 
which is simply proportional to the dynamical timescale in the 
disk, too many galaxies are formed, and the number of galax- 
ies at faint-end is greatly overpredicted. We found that the SF 
timescale should be nearly constant against redshift, as sug- 
gested by our previous analysis (NTGY) and by other recent 
SAM analyses(Kauffmann & Haehnelt 2000; Somerville, Pri- 
mack & Faber 2001). 

The uncertainty in estimation of galaxy size results in an 
uncertainty in estimation of the surface brightnesses of galax- 
ies, which is directly related to the selection effects mentioned 
above. In our SAM, like usual SAMs, the size of the disk is de- 
termined under the assumption of specific angular momentum 
conservation of cooling gas, so that the disk size is proportional 
to the cooling radius. In order to see how theoretical predictions 
are changed by changing the galaxy size, we introduced a free 
parameter p allowing for an additional redshift dependence in 
size estimation. We found that the value of p > 0.5 is favored 
in order to reproduce the observed counts, which indicates that 
the disk radius should be extended by a factor of (1 + z) over 
the cooling radius. We also found that this manipulation cannot 
be discriminated observationally because we cannot know how 
many undetected low surface brightness galaxies below the se- 
lection criteria are there at high redshift, which are presumably 
systems of large size. 
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Through this work, we have shown that our SAM can explain 
a variety of observed properties of nearby and high-redshift 
galaxies, and place some constraints on star formation, size 
evolution, and dust extinction. More stringent constraints will 
certainly be obtained by a greater knowledge of dynamical and 
kinematical properties of galaxies in near future. 
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Fig. 5. — Number-magnitude relations for various cosmological models in the Subai'u K' -hand. The solid, dot-dashed, and dashed lines indicate SC, OC, and LC, 
respectively. The thick lines denote the models including the selection effects, while the thin hues denote those without the selection effects. The open circles with 
etTorbars indicate the raw counts in the SDF. The other symbols indicate other observational data refen'ed in the figure. 
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Fig. 6. — Number-magnitude relations for various cosmological models in the HST U BV I-bwds. Types of lines are the same as Figure 5. The symbols indicate 
observational data referred to in the figure. 
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Fig. 7. — Isophotal area of galaxies detected in both the K' and / bands in the LC model. The parameter p is varied from (solid line), via 0.5 (dashed line), to 1 
(dot-dashed line). The la scatters in the model with p = are shown by error bars. The crosses indicate the data in the SDF. 
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Fig. 8. — Redshift distributions of galaxies for a A-dominated flat CDM model. Left panels: the /814-band. Right panels: the ^T'-band. The soHd and dashed lines 
denote the LC model with and without the selection effects, respectively. The histograms in the left panels indicate the redshift distributions of the HDF galaxies 
derived by Furusawa et al. (2000), using an improved technique of photometric redshift estimation. 
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Fig. 9. — Number-magnitude relations for a A-dominated flat CDM model. The solid, dashed, and dot-dashed lines indicate the LC (CSF), LDl (DSF with a = 1), 
and LDO (DSF with cr = 1) models, respectively. 
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Fig. 10. — Number-magnitude relations for a A-dominated flat CDM model. The solid, dashed, and dot-dashed lines indicate the LC models with different 
z-dependences of galaxy size, p = 0, 0.5, and 1, respectively. 
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Fig. 11. — Number-magnitude relations for a A-dominated flat CDM model. The solid, dashed, and dot-dashed lines indicate the LC models with different 
2-dependences of dust optical depth, 7 = 0, 1 , and 2, respectively, of the slab dust model. 
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Fig. 12. — Number-magnitude relations for a A-dominated flat CDM model. The solid, dashed, and dot-dashed lines indicate the LC models with the slab dust 
with 7 = 0, and the screen dust with 7 = and 2, respectively. 




15 20 25 20 25 30 



K' magnitude . magnitude 



Fig. 13. — Morphology-dependent number-magnitude relations in the LC model. The solid and dashed lines indicate the total and early-type (E/SO) galaxy counts, 
respectively, in the LC model. The thick lines of these two are for /bulge = 0. 1 , and the thin lines for /bulge =0.5 model. 
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Table 1 



Model Parameters 



CDM Model 


co^ 


;mological 


1 parameters 




astrophysical parameters 






Ho 




h 




Vhot(kms ') 


Ohot 


r", (Gyr) 


a* 


A 


T 


SC 


1 





0.5 


0.67 


350 


5 


7 


-3 


2 


1 


oc 


0.3 





0.6 


1 


240 


3 


L8 


-2.2 


1 


2 


LC 


0.3 


0.7 


0.7 


1 


280 


2.5 


2 


-2 


1 


L5 



Note. — Equations defined the SN feedback-related parameters (Vkot and ahot) and the SFR-related parameters (t° 
and a») are (2) and (1), respectively. The definitions of last two parameters, fi, and T, are given in §§2.3 and 2.2, 
respectively. The other parameters used in common as standard include Qb = O.OlSh'^, CSF (equation 1), /Mge = 0.1 
(§2.1), 7 = (equation 3), p = (equation 6), and (A, (tx)=(0.05, 0.5) (equation 5). 



